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Abstract. We discuss two recent applications of homotopy analysis method and 



homotopy perturbation method and conclude that the results are completely useless 



from both mathematical and physical points of view. 



In two papers appeared recently in this journal El-Wakil and Abdou [1] and 
Abdou [2] proposed the application of homotopy analysis method (HAM) and homotopy 
perturbation method (HPM), respectively, to several nonlinear equations of supposed 
interest in fluid mechanics, plasma physics, optical fibres, biology, solid state physics, 
chemical kinematics, chemical physics, geochemistry [1], particle vibration in lattices, 
currents in electrical networks, pulses in biological chains, etc. [2]. 

These two papers are examples of the new physics and mathematical physics that 
has lately spread among many journals. In what follows we discuss them. 

El-Wakil and Abdou [1] stated that "Being different from the perturbation 
technique, the HAM does not need any small parameter to function". However, if 
we look at the paper we realize that they expanded a solution to the nonlinear equation 
in the way that is common to perturbation theory: u(x, t) = u (x, t) + Ui(x, t)p + . . . + 
u n (x, t)p n + . . . and set the dummy perturbation parameter p equal to unity at the end 
of the calculation. What is p if not a small parameter?. Most curiously, in the other 
paper Abdou [2] stated that "Using the homotopy technique in topology, a homotopy 
is constructed with an embedding parameter p e (0, 1) which is considered as a small 
parameter" . It seems that a dummy perturbation parameter that plays exactly the same 
role in both approaches is a small parameter in HPM [2] but it is not in HAM [1]. 

El-Wakil and Abdou [1] studied the following system of coupled nonlinear 
equations: 



du 

~dt 

dv 

dt 



u(l — u — v) + 
d 2 v 



d 2 u 



dx 2 




(1) 



with the initial conditions 



—kx 



u(x, 0) 



(1 + e ~ fc */ 2 ) 2 



HAM and HPM: Another attack to reason 3 

^°) = 1 + e L/ 2 (2) 
where A; is a constant. According to the authors "The above system of two coupled 
nonlinear equations of reaction diffusion type arise in chemical reactions or in ecology, 
and other fields of physics". It is remarkable that they are able to solve problems in 
so many fields of scientific interest with such simple ubiquitous system of equations. 
Unfortunately, the authors did not give any reference about such applications. 

They applied HAM and derived a most cumbersome approximate solution that they 
displayed in their Appendix A. If you look at it you realize that it is merely the set of 
Taylor expansions 

u(x, t) — 22 u j( x W 

3=0 

v(x,t) = ^^Vj(x)t 3 (3) 

3=0 

through third order. Obviously, if you substitute these equations into equation flTJ you 
easily obtain the coefficients Uj(x) and Vj(x) that the authors so laboriously derived by 
means of HAM. 

According to El-Wakil and Abdou [1] the exact solutions are 

U ^ = n I kz/2\* 

(1 + e Kz i l ) 

where z = x + xt. However, the reader may easily verify that they do not satisfy 
the nonlinear equations (CQ). The authors stated that those exact solutions come from 
another paper [3] where one reads 

U(z) = ; — - 

{l + e kz/2f 

and z = x + ct. Unfortunately, the second equation does not satisfy the corresponding 
initial condition and, therefore, it cannot be a solution. It seems that the system of 
nonlinear equations <^ exhibits an exact solution but the authors failed to provide it 
in two papers were they treated the same problem with two different approaches [1,3]. 
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Abdou [1] discussed the differential difference equations (DDEs) 

du r 



(l + au n + (3u 2 ^j (u n+ i - u n -i) (6) 



Ot 

where a and (3 are constants. The author did not say it explicitly but it appears that 
n is an integer. According the author "The hybrid nonlinear difference equation ([6]) 
describes the discretization of the Korteg-de Vries (KdV) and modified KdV equations 
with initial condition" 

Un,o — a o {aa + 2) tanh 2 (A;) tanh 2 (/cn + c) (7) 

"whose exact solution reads as" 

1 



u n {t) = ao (aao + 2) tanh 2 (A;) tanh 2 

a 

x kn + (aao + 2) 2 tanh(A;)sech 2 (fc)- + c 



(8) 



They applied the HPM and obtained just the first term of the Taylor expansion of u n (t) 
about t = 0. The author also treated other DDEs with known exact solutions. Since 
they are simpler than the one discussed above he obtained some more terms of the 
Taylor series. 

Summarizing: El-Wakil and Abdou [1] and Abdou [2] applied HAM and HPM, 
respectively, and obtained Taylor series of differential equations with known exact 
solutions. Such equations are merely tailor-made toy problems with the purpose of 
applying those methods. One easily obtains the time series by straightforward expansion 
of the equations about t = 0. The resulting approximate solutions are valid for 
short times, which are typically of little or no physical interest at all. Although the 
authors mentioned that the models appear in every branch of physics, chemistry, and 
engineering, it is clear that they do not apply to actual problems of interest. 

Variational and perturbation approaches (VAPA) like HPM, HAM, Adomian 
decomposition method (ADM), variational iteration method (VIM), etc. have produced 
many pseudo scientific works of the poorest quality. We have denounced this unhappy 
situation in several articles [4-16] but apparently nobody seems to care. We suggest 
the reader to have a look at some of the articles cited by El-Wakil and Abdou [1] and 



HAM and HPM: Another attack to reason 5 

Abdou [2] in order to have an idea of what we are talking about. There is a sort of 

endogamic refereeing process that spreads VAPA like vermin and more journals are 

being contaminated every day. 
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